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In Yang-Mills theories, variational calculations of the Rayleigh-Ritz type face the problem that 
on the one hand, calculability puts severe constraints on the space of test wave functionals; on the 
other hand, the test wave functionals have to be gauge invariant. The conflict between the two 
requirements can be resolved by introducing a projector. In this paper we present an approach 
to approximating the projector in a way known and successfully employed in nuclear physics: the 
Kamlah expansion. We discuss it both for electrodynamics and for Yang-Mills theories to leading 
order in a perturbative expansion, and demonstrate that the results are compatible with what one 
would expect from perturbation theory. 
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INTRODUCTION AND MOTIVATION 



Recently, there has been a renewed interest in the application of variational techniques to Yang-Mills theories, cf. 
I [Q, H, ^ 1^, 1^ and references therein. Whereas the first investigation in this line Q did not focus explicitly on the 
subject of gauge invariance, subsequent investigations have put more stress upon the subject. The basic idea can 
^\ [ already been found in the nuclear physics environment if one starts from a trial wave functional that does not 
share the basic symmetries of the Hamiltonian, one should not vary the expectation value of the energy but rather 
[ introduce a projection operator P that projects the wave functional onto the symmetric space and vary this new 
. expectation value: 

0-( This construction allows the usage of simple wave functionals. It turns out, unfortunately, that the projector acting 
on e.g. the simple class of Gaussian wave functionals cannot be evaluated exactly in the case of Yang-Mills theories. 
^ Therefore a number of approximation methods have been devised. One class of approximations treats the projection 
■ integral as an effective non-linear sigma model (with a non-local action) and uses common field theoretical approx- 
^ ■ imation methods to deal with this model ||], cf. also [||. Another approach has suggested an 'improved energy 
•'H functional' jj, ||. The idea that is followed in this approach is that the presence of the projector could be simulated 
' by considering an improved energy functional which contains besides the energy expectation value additional terms 
^ , which depend on the generator of the symmetry (x) : 

E-ff = {i;\H\^) - {e-^r\K,y){^\T-{^)r\ym (2) 

where we have assumed that the states are properly normalized. Here O"* denotes the "moment of inertia" corre- 
sponding to the "rotations" generated by r"(x). In this paper we want to investigate how one can derive such an 
improved energy functional and under which restrictions this can actually be done. 

The structure of this paper is as follows: first, in section ||, we present the projector onto gauge invariant states as 
it is be used in this paper; we also comment shortly on how one can write down such a projector in the presence 



of external colour charges. In section III we introduce the general concept of the Kamlah expansion as it is used in 
nuclear physics 0, ^ , and formulate it in a way s.t. it may be used also for gauge theories. In section ^ we 
apply the concept to the theory of a quantized electromagnetic field in the presence of classical charges. There we 
will see that the leading order of the expansion - which was often sufficient in the case of nuclear physics - is in fact 
inconsistent in quantum field theories. In section|^we discuss the application to Yang-Mills theories. However, we are 
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only able to discuss the leading order of a subsequent perturbative expansion. Also, some additional difficulties w.r.t. 
the usage of external charges will become apparent. In section ^ we give a summary and a critical evaluation of the 
Kamlah expansion. Conventions and the canonical quantization scheme used in this paper are discussed in appendix 
some results from perturbation theory will be presented in appendix ^ and details of some of the calculations are 
given in appendix 

II. THE PROJECTOR ONTO GAUGE INVARIANT STATES 

A. Generalities 

In this section we want to introduce the concept of the projector onto gauge invariant states. The basic idea is 
simple: if we start from a wave functional ■i/;ni[A] = (AjV'ni) which is not gauge invariant, we can construct a gauge 
invariant wave functional ^/'gi[A] by summing over all gauge transformed wave functionals with the correct Haar 
measure. Symbolically 



^gi[A] = / PA*[0]^ni[A^] 



(3) 



symbolizes the integration over all (as we will see below: topologically trivial) gauge transformations 
For notation cf. also app. |A[ For the purposes of this paper - as outlined in app. A4 



it IS 



where 

parametrized by 

sufficient to have wave functionals that are annihilated by the Gauss law operator. This is equivalent to having wave 
functionals that are invariant under topologically trivial gauge transformations; therefore the domain of integration in 
eq. (§) can be restricted to topologically trivial gauge transformations also. In this case the gauge transformation can 
be written with the help of the Gauss law operator Q (in ab senc e of external charges given in eq. (^, in the presence 
of external charges a further discussion can be found in sec. II B ): 



V'ni[A] 



This now allows to rewrite eq. (H) in terms of an operator PP acting upon ipn 



^gi[A] = / P^[<^] 



-ij (i^a;0''(x)e"(x) 



= P^VniiA], 



(4) 



(5) 



where denotes the integral over all variables parametrizing the gauge transformation including the Haar measure 

necessary to make the integral invariant under left multiplication. 



B. Projector in the presence of external charges 



If we restrict ourselves to the pure gauge sector without external charges the Gauss law operator is simply (for 
conventions and placement of factors of g cf. appendix is defined in eq. ( AlCl| )) 

g'^(x) = -F'^(x) (6) 

Both the Gauss law operator and its pure gauge part fulfil the "angular momentum" algebra 

r (x), g\y)] = zr^^g^(x)J.y, [-F'^(x), -F''(y)] = zr''=(~F^(x))5.y . (7) 

If we want also to consider external charge contributions, this can be done most economically in the following way (fo r 
definiteness, we assume the external charges to be fermions, in the following gcnerically called 'quarks') ||, [ll[ : 
add to the Hamiltonian of the gauge sector a matter field part Hm- 

(8) 



Hm = M J d^xV'(x)?A(x) 

where AI denotes the quark mass, and we treat a single quark species ("flavour") for simplicity. In concordance with 
the common interpretation we associate creation operators with operators that create positive energy states, and thus 
write explicitly 



C(x) = 



/ a?(x) \ 
ai(x) 
6^(x) 



V'r(x) 



4''(x) 

&?(x) 
V &f(x) / 



(9) 
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We then interpret the operators a, as connected to quarks, and the operators 6, as connected to anti-quarks. We 
now turn to the charge operator. Let A" be the (hermitian) generators of the SU(N) representation appropriate for 
the quarks. The charge operator is then given by 

p"(x) - Jo^x) = ^(x)70AXx), (10) 

where there is no integration over x just as in the remainder of this section (unless exphcitly stated). In terms of the 
a, b operators, it can be written as 

p"(x) = al''(x)A^,<(x) - 6l''(x)A^:6^(x). (11) 

The matrix A* denotes simply the complex conjugate matrix of A. The charge operators inherit their commutation 
relations from the A matrices: 

[p''(x),p''(y)] = 5xy (al(x)[A^A'']a.(x) +6l(x)[A-,A''1&.(x)) = zrV(x)<5xy. (12) 

Thus, we can form the Gauss law operator in the presence of external charges 

e'^(x) = -r''(x)+p''(x), (13) 



where r"(x) is defined in eq. (AlO). We form states in the presence of external charges (called jcoupled)) in the 
following way: 

Icoupled) = ^ c°|fermion i)|YM a), (14) 

where are probability amplitudes and |YM a) are purely gluonic states labelled by a. The fermionic states are 
generated by acting with quark and anti-quark creation operators on the (gauge invariant) fermionic vacuum ^: 

Ifermion z) = al ■ ■ ■ b]^ ■ ■ ■ b]^ |0) , (15) 

where we have combined colour-, spinor- and spatial indices into a single super- index i (or j), and denoted the 
fermionic vacuum by |0). We can also write suggestively 



a,ii:---,3m. 



•«l4---^L|0>|YMa) = «l---«l4---dO)|YM),,...,„ (16) 



with |YM),,...,„ = Ea<i...,,„|YM a). Using 

e-^i ''"y^''(y)p"(y)i:''^{^)e'S ''"y^"^y)p"(y) = (pVW^'')„cC(x) - qf,[C/(x)]C(x), 
e-^S d,'yv\y)p''(y)^l^{-i^)e-! d'yv\y)p''{y) ^ ^tc(x)(e-»'^''W^''),„ = pS;[t^(x)]^i"(x) 

where transforms according to the representation q^^[C/(x)] and t/;^ transforms according to the conjugate repre- 
sentation p^^[J7(x)] — q^[L/^^(x)], the requirement 



coupled) — Icoupled) (18) 



leads to (basically) the same result as known in the literature where one considers the Yang-Mills wave functional 
|YM)ij.. with the extra colour indices: 

g./.3.0»(x)r"(x)|YM),^...^.^^ =q« jC/]...q«^J[/]p«^Jf/]...pj?^^^[[/]|Y^^^ (19) 

where the super-index notation for the representation matrices q,p has to be understood as follows: if we decompose 
the super-index ii into the colour index ai, the spatial index mi and the position index xi, and correspondingly li 
into 6i, ni and xi (the position index is in both cases identical since we do not associate an integration over position 
indices with gauge transformations), then 

{q^,^%,iAU] - <5™,„,(q^,p«),,bjC/(xi)]. (20) 

The formulation chosen in this section has the advantage that one can give a closed expression for the projector even 
in the presence of external charges. At this point one has to note that from the fact that Q°'\ip) — one may not 
conclude that one has diagonalized all generators of the gauge part of the Gauss law operator, since the matter part 
pa ga ^^sQ Operator valued. One can only choose the Cartan subalgebra of the colour charges to be diagonalized. 



^ That the fermionic vacuum is gauge invariant is a trivial consequence of the fact that it is annihilated by the annihilation operators. 
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C. The Projector onto Physical States 

Since we now have an operator that generates gauge transformations on arbitrarily charged states, we can also 
write down a very simple (formal) expression for the projector onto the physical sector: 



The integration domain runs over the topologically trivial sector of the gauge transformations. Two points deserve 
attention: first, we want to consider why PP is indeed a projector onto the physical subspace, and second, why it is 
sufficient to restrict the domain of integration to topologically trivial gauge transformations. In fact these two points 
are closely related. A state in the physical subspace is characterized by the fact that it is annihilated by the Gauss 
law operator: Q'^l'ip) = 0. An identical requirement is 'f'"^" {ip) — lip), in other words a state has to be invariant 
under topologically trivial gauge transformations. We have abbreviated here / cf x (f>°' {x)Q°' (x.) as <j>'^Q"'. If we now 
act with e"' '^ ^ upon P'' we can use (a) the fact that topologically trivial gauge transformations form a group, i.e. 

(b) the left invariance of the Haar measure 

Vfi[ip]^Vfi[a{cj),^)], (23) 

and (c) the fact that a topologically trivial gauge transformation does not change the domain of integration of the 
projector P'' to obtain 

= J Vfi[ip] e-'J "'"W^'/')^" ^ J x>ti[a] e-'^""^" = P". (24) 

Thus we have answered both questions: the projector projects onto the physical subspace since the projector is 
invariant under small gauge transformations, and since the Gauss law constraint is equivalent to topologically trivial 
gauge transformation only this domain of integration is required. 

III. GENERAL CONCEPT OF THE KAMLAH EXPANSION 

A. Abehan Case 

The concept of the Kamlah expansion Q is most easily explained using an Abelian quantum mechanical example 
1^, 0, 1^; assume that we have a two-dimensional system. The generator of rotations is J. The projector onto states 
having angular momentum / is 

pi = r ^e^'^iJ-n, (25) 







27r 

That this is indeed the correct projector can again be seen most easily by 



g»/3./p/ ^ '^^i(a+0)J-i(c+l3)I+ipi ^^ipipl /26) 

io 2^ 

where it has been used that = J^^^^ due to the periodicity of the integrand. 

The case where the Kamlah expansion works best is the case of strong deformations. Strongly deformed states 
are states where 

Me^"-'m (27) 
decreases rapidly for increasing a. This is very useful, since one can argue that 

MHe-^'m (28) 
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will vanish rapidly for increasing a, too, but 



(29) 

is a smooth function of a. The argument consists basically of noting that, whereas J is a one-particle operator (in 
many-body language), e*""^ is a collective operator that affects arbitrarily large numbers of particles 

6"-^= ^ + iaJ + ^{iaJ)^ + ^{iaJ? +■■■, (30) 

O—body operator 1 — body operator ^ v ^ v 

2 — body operator 3~body operator 

H is (in nuclear physics) at most a two-body operator (A similar terminology can also be constructed for Yang-Mills 
theories, which is done in In this terminology the Yang-Mills Hamiltonian also will be at most a two-body 

operator, i.e. contain up to four creation/annihilation operators.); thus the behaviour w.r.t. increasing a should be 
qualitatively the same whether we consider e'"'' or He'"''. One can now argue further that in the matrix element 
{(f)\He'^°''' \(t>) all degrees of freedom that are not affected by the collective rotation ('internal degrees of freedom') 
can be integrated out, and that this results in an effective Hamiltonian operator that is an - up to now - unknown 
function of J, the symmetry generator. For this function it is proposed to use a power series expansion in powers of 
(J — (J)) (where (J) is used as an abbreviation for (0| J|(/))), which is plausible, since one does not expect singularities 
for collective rotations with finite angular momenta. (The idea to expand the Hamiltonian in powers of the symmetry 
generators, although with a different method to obtain the coefficients, was proposed in cf. also ^.) 

However, this is only useful if we can stop the expansion after a few terms. That this is indeed possible one can see 
as follows: 

• by assumption, (0|e*"'^|0) is well localized in a space, since we only consider strongly deformed states; therefore, 
a broad range of wave functions with good J quantum number have to be added up coherently 

• ((/>|( J— (J))e*"'^|0): each wave function with good J quantum number is now weighted by that quantum number 
and this begins to destroy the coherence, in other words (0|(J— (J})e'"'^|(/)) is broader in a space than (0|e*"'^|(/)) 

• the higher the power of (J — {J)), the broader is the resulting matrix element in a space, hut we already know 
that {(j)\He^°'-'\(j)) is not very much broader in a space than (0|e'"''|</)), hence the matrix elements containing 

higher powers of (J — (J)) (besides e*"'') must be suppressed (We assume here that there is no 'conspiracy' 
between different orders in the expansion, but this seems to be a plausible assumption.), i.e. the corresponding 
expansion coefficients must be small, qed. 

Thus, we can stop the expansion 

(<^|i/e-'''|0) = VA,(0|(J- (J))'^e»'^"|</.) (31) 



after a few terms: 



k=0 



\Happro.e'"'\<l>) = ^Afe(0|(J- (i))'=e^"^|0). (32) 

fe=0 



The (n-|-l) coefficients are determined by requiring that {(j)\H e'"'^ and {(t>\Happroxe-'^°''^\(t>) along with their first n 
derivatives w.r.t. a are equal at the point a = 0. One thus ends up with the set of equations^ 

m{j-{j)n) - ^ A,(</-i(j- (j))'=+™i0), (33) 

fe=0 



^ In order to obtain a more compact notation, we have equated (0|He"*'''~^''^' |0) and {<j)\Happroxe^°'^''~^''''^\4>) ■ In the sections in 
field theory, we will equate the terms equivalent to (0|lfe"^(''~^' |0) and {4>\Hapj}Tox^^°'^'''~^^\<l>)- They difl'er only by a multiplicative 
c-number e""^^"'^^ or e~"^^ respectively. 
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where m = 0, . . . , n. After we have solved these equations to obtain Ai, we can insert Happrox = Sfe=o ^k{J ~ {J))'^ 
into the expression for the projected energy 

The advantage of the approach now becomes visible: whereas before we did not know how to evaluate the matrix 
elements including the projector, it has now become trivial; P-^ projects onto angular momentum I. Then JP^ = IP^ , 
and numerator and denominator cancel, which means that we do not have to calculate a single matrix element 
containing the projector explicitly: 

n 

E'a7pLM) = T.Mi-{jyf- (35) 

k=Q 

In nuclear physics, for most applications n is taken to be one or two; one then ends up with 

E^p-pioM) - + jf^il (J)) (36) 

{{Ajy} 

for n = I and 

Erp-pioAi) = {H) - + ^(/ - (J)) + - (J))^ (37) 

with J- = ^^'^-^r^ and^^i:^ ^ ^ (38) 

2Iy (AJ4)-(AJ2)2 Isc (AJ2) 

for n = 2 We have abbreviated (0| . . . by (. . .) and used the notation AJ — J — (J). If we take simply n — 0, 
we end up with Aq — {H) which is reassuring: the lowest order term just gives the unprojected expectation value. 
We see that, at first order, we get no corrections to the mean field energy if (J) = /, whereas for the expansion to 
second order, we even then obtain a correction. One can interpret this correction if one considers a deformed state as 
being a superposition of states with different angular momenta; especially, there are components in the wave function 
that are of higher angular momentum than the angular momentum one projects onto. Usually one would assume that 
higher angular momenta also contain higher kinetic energy that should not be present in the projected wave function; 
thus one obtains a negative correction to the ordinary energy expectation value. One notes in passing that the two 
moments of inertia Iyj^sc are quite different from one another; the latter is under certain circumstances identical 
to the moment of inertia of the cranking model Q, since one can show that, if the variational space considered for 
the cranking calculation contains not only |0) but also AJ\(f>), then the solution of the cranking equations indeed 
solves also the variational equations derived for the n = I Kamlah expansion, and that the Lagrange multiplier of the 
cranking model has precisely the same form as the first order correction of the Kamlah expansion p9| 
There are calculations that implicitly assume that in practice the moment of inertia ly can be determined by a 
cranking type calculation ||, |^. In Yang-Mills theories it turns out that to O (<?'') this is indeed the case, but the 
general expressions given above in eq. (|38|), and the fact that the variational space of Gaussian wave functionals is to 
small to prove the equivalence of Xsc and the cranking model moment of inertia, do not give much hope that it stays 
this way for higher orders in perturbation theory. 



B. Non-Abelian Case 



If one goes from an Abelian type problem to a non-Abelian problem (like three dimensional angular momentum 
projection) one immediately has to face some problems ||2^; the first one is that a rigid body has two frames of 
reference (the lab frame and the intrinsic frame), two sets of angular momentum operators and (since the total 



^ The variational space that we use in field theory, consisting of Gaussian wave functionals only, is not so large as it is required for this 
equivalence. 
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angular momentum is identical in the intrinsic and in the lab frame) three angular momentum quantum numbers 
(usually called I,K,M belonging to — P , J3, I^) The basic ingredient into any projector is then 

P'km = \I,K){I,M\. (39) 

The quantum numbers I,K are those of the state to be observed in the lab frame, whereas M is only seen in the 
intrinsic frame. The object Pj^j;,/ does not have the properties usually expected of a projector, since it is neither 
hermitian, nor do we have a relation like = P. The obvious remedy cannot be used, since taking P^^^j ^.^ would 
avoid the mathematical problems but would violate rotational symmetry, as an 'intrinsic' state with M — Q would 
always lead to a state to be observed in the lab frame with X = |P . One could also ask more physically: if M is an 
internal quantum number that has rather little to do with what can be observed in the lab frame, why don't we use 



Pk = Y.f^tPKM (40) 

M 



as a projector with Jm as additional variational parameters (cf. e.g. |^) ? In it was argued that one should take 
all /a/ to be equal to one, but whatever one does, one always ends up with the problem that one has to evaluate 
matrix elements of operators acting upon projectors that do not carry the quantum numbers that are projected upon; 
to be concrete, one has e.g. to evaluate 



Em.m'{^\PL 



(41) 



MM' 



One sees clearly that, unlike in the Abelian case, the projector matrix elements do not usually cancel. A number 
of ways have been devised in the literature to approximately evaluate these sums, but they are not of interest here, 
since they cannot be transcribed to the field theory case. We only want to point to the fact that the difficulties in 
even the formulation of the projector do not come into play if we want to project onto a state with / = 0, since 
then there is only one projector Pqq and no factors of / can be chosen. This is the case of the Yang-Mills projector 
since we have combined the gluonic (F) and the matter (p) contribution into one set of operators (Q) which now are 
supposed to form an 's-wave'. The problem of non-factorizability of projectors and other operators will nevertheless 
also haunt us later on, too (cf. sec. |V A ). 



C. General Remarks about Validity 



One last comment shall be made at this point on the validity of the Kamlah expansion. It is a priori not clear 
that it is physically sensible to perform the Kamlah expansion for a local symmetry. In nuclear physics, the very 
deformed states usually do not come about through a very deformed nucleus in space, but are due to a very large 
number of participating nucleons In a gauge field theory, one has usually the situation of a very small number 
of excitations at every point in space. Therefore it is not clear whether it is sensible to use the 'large-deformation' 
expansion. However, it will turn out that the Kamlah expansion performs very well in the cases we have considered. 



IV. APPLICATION TO QUANTIZED ELECTRODYNAMICS 

In this section, we want to apply the general principles outlined above to what we call electrodynamics - a quantized 
electromagnetic field interacting with static sources. It is illuminating since we will see that in marked contrast to 
nuclear physics, the first and second order expansions are even qualitatively quite different. 



A. Space of Allowed Wave Functionals 

In the remainder of this paper we will restrict ourselves to wave functionals of Gaussian form 

^P[A] = exp |-(A - A), Qg,--1 - iEy^ (A - A), + ze,(A - A),| (42) 

where we use a super-index notation: when we consider electrodynamics, i,j both include a spatial and a position 
index, which is summed or integrated over, respectively. In the case of Yang-Mills theories the super-indices further- 
more include a colour index. 0"^,^, A,e are all purely real. Further conventions and a short account of canonical 
quantization are given in appendix H. 
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B. Simplified Projector 

In the case of electrodynamics the expression for the projector can be simphfied significantly compared to eq.(pTl). 
Since in electrodynamics there is only one charge operator, which also commutes with the Hamiltonian, the fermionic 
states can be chosen as eigenstates of the charge operator, and since we want to prescribe a certain external charge 
configuration we consider only states (cf. eq. (|4|) for the non-Abelian case) of the form 

|coupled) — jfermion p}|photon), (43) 

with 

/3(x)|fermion p) = /5(x)|fcrmion p). (44) 
Since we only want to project operators that depend solely on the gauge field part, we can use a simplified projector 

P''|coupled) y"p0e'^'^(^-'')|fermion p)|photon) = Ifermion p) ^ P0e'^ '^(^"'^^Iphoton) 

= jfermion p}P''| photon) (45) 

which does not depend on the fermion charge operator but directly on the prescribed c-number external charge 

distribution. Furthermore, the (pure gauge part of the) Gauss law operator is also simpler than in Yang-Mills 
theories. It is given by 

r(x) = v,n,(x). (46) 



C. Kamlah Expansion to First Order 

We start from the expression (cf. eqs. (^, |3^)) 

(iJe'/0(r-P)) ^ Ao(e'^*(r-''))+Ai(y)({r(y) - (r(y))} e'/^^^-p)^^ 



(47) 



where we have again used J (j)(r — p) as an abbreviation for J (Pz 4>{z){T{z) — p{z)). A useful notation that will be 
used throughout the remainder of this section is 



A(x) = r(x) - (r(x)) , 5(x) = p(x) - (r(x)). 

Introducing the functional inverse 0^^(x, y) of (A(x)A(y)}, s.t. 

e-i(x,z)(A(z)A(y)) =5xy 



(48) 



(49) 



we obtain for the projected energy functional (details of the procedure in general are outlined in app. where it is 
applied to the specific case of the second order Kamlah expansion) 



(HPP) 



(PP) 



(i7) + (i7A(x))e-i(x,y)5(y), 



(50) 



where we have introduced the notation Ajj^at to denote the fact that A is evaluated using the Kamlah expansion to 
first order. Note that the derivation can be carried out for any gauge-invariant few-body operator O, since there is 
nothing special about the Hamiltonian: 



(OPP) 



{PP) 



= (0) + (0A(x))e-i(x,y)5(y). 



(51) 



This provides a first test of the validity of the first-order Kamlah expansion. We may set O ~ A(z), and obtain 



(A(z)PP} 



I.e. 



(PP) 

(r(z)P^) 
(pp) 



eq. (51) 



■4 (A(z)) + (A(z)A(x)} e-i(x,y)<5(y) = p(z) - (r(z)). 



eq. (til 



(52) 
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showing that the expression provided by the first order Kamlah expansion for (r(z)P'')/(P'') is exact. It doesn't stay 
that way, however, if one considers e.g. O = A(zi)A(z2) : 



(A(zi)A(z2)P'') 



eg. {51) 
1st 



A(zi)A(z2)) + (A(zi)A(z2)A(x))e-^(x,y)<5(y) - (A(zi)A(z2)), (53) 



where we have used that (A(zi)A(z2)A(x)) — 0, which is true for every Gaussian state as given by eq. (^2|). We see 
that there are no corrections due to the projector for A(zi)A(z2) by the first order Kamlah formula. This gives rise 
to the expectation that the energy will not be correctly projected, since it contains terms of the type (r(x)r(y)P'') 
(as we will see below, they arise from the kinetic energy). 

Let us now evaluate the projected energy. For this it is useful to note that under a variety of assumptions'* 

(Bi(x)B,(y)r(z)) = (B,(x)B,(y))(r(z)), (54) 

i.e. (Bi(x)Bj (y)A(z)) = 0. Physically it is plausible that the magnetic part of the energy should not be affected, 
since it contains the transversal degrees of freedom only. For the electrical field one obtains 

(n,(x)n,(x)r(z)) = 2(n,(x)A(z))(n,(x)) + (n,(x)n,(x))(r(z)}, (55) 

which is valid even without integrating over x. With this, the first order energy correction is given by 

{HA{z)) = (n,(x)A(z))(n,(x)). (56) 
One can make further progress by assuming that both G^^ and S are translationally invariant; then one can write 

(n,(x)n,(y)), = iG-.i(x,y) +4(EGS),,-(x,y) = J ^e^P-(-y)(aPpf/ + a^P^--) (57) 

which defines the longitudinal and transversal components fi^, ct^, and where we have used the abbreviation Pf^^ for 
PiPi/P^; correspondingly P^" — 5ij — Pf ". This allows an exphcit expression for 9^^(x,y): 



(58) 



and thus 



(i/A(xi))e-i(xi,y) = (n,(x))(n,(x)A(xi))e-i(xi,y) 

= (r(x)}GA(x,y). (59) 

Ga denotes the Green's function of the (ordinary) Laplace operator, i.e. — ViViGA(x, y) — (5xy. Then the projected 
energy reads 



{HPP) 



(PP) 



^= {H) (r(x))GA(x,y)((r(y)) -p(y)). 



(60) 



We can make this expression even more transparent if we have a look at (Ilf (x))(nf'(x)) where Il f (x) denotes the 
longitudinal component of ni(x), i.e. nf(x) = / d^y {P^)ij{x,y)Ilj{y); for (P-^)ij(x,y) cf. eq. ( A19 ). One can 
derive the identity 

(nf (x))(nf (x)) = (r(x))GA(x,y)(r(y)), (61) 

and therefore one can write 



{HPP) 



{PP) 



{Hf - i{(r(x))GA(x,y)(r(y)) -2(r(x))GA(x,y)p(y)} (62) 
{Hf ~ \ (r(x) - p(x)) GA(x,y) (r(y) - p(y)) + ip(x)GA(x, y)p(y). (63) 



* Either one may assume that the kernel in the Gaussian wave functional is purely real, i.e. S = 0, or that both G ^ and S arc translation 
invariant, or that the background field = EijTjVjAfe fulfils the classical equations of motion. 
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{H)^ denotes here the part of the energy expectation value which is independent of (H^), but which still depends on 
(n^n-^)c, where we have introduced the abbreviation {AB)^ = {AB) — {A){B). In formulas this is expressed simply 
by (H)^ = (H) — i(nf (x))(nf (x)). Thus (r(x)} is determined to be equal to p(x) at the stationary point, giving 
the correct charge contribution to the energy 

ip(x)GA(x,y)p(y). (64) 

As can be seen by explicit computation, at the stationary point (n^n^)c is zero, s.t. the energy depends only on 
the transversal degrees of freedom and on the charge density. However, the outcome is quite in contrast to what we 
expected from the Kamlah expansion initially: since we have computed the projected energy functional, (a) this should 
be completely independent of the gauge dependent degrees of freedom; thus, it should not depend on (n^n-^)c, since 
in order to provide a non- vanishing expectation value (n^n^)c the wave functional has to depend on A^; and (b) 
the state was assumed to be strongly deformed, but, if we take the first-order formalism seriously, at the stationary 
point our state is just not deformed, since (n^n^}c = 0. These problems are overcome by performing the Kamlah 
expansion to second order. Therefore, we now turn to this. 



D. Kamlah Expansion to Second Order 



properties of Gaussian wave functionals, cf. eq, 
in this derivation. We obtain for the projected energy functional (using the notation AH = H — {H)) 



Details of the derivation of the energy functional are given in app. CI, It has to be noted that the factorization 

and translational invariance of both G^^ and S are used heavily 



(65) 



{HPP) 



(PP) 



2nd 



{H) - -(AifA(xi)A(x2))e-i(xi,x2) 

-(AifA(xi))e-i(xi,yi)(r(yi) - p{y,)) 

+ l(Ai/A(xi)A(x2))e-i(xi,yi)e-i(x2,y2)(r(yi)-p(yi))(r(y2)-p(y2)), 

where we have introduced the notation A\r^„a to denote the fact that A is evaluated using the Kamlah expansion to 
second order. As in the paragraph above, in the derivation of eq. ( |65| ) no special property of H apart from gauge 
invariance was used, and therefore the expression is valid for every gauge-invariant few-body operator O: 



(OPP) 



(PP) 



2nd 



{O) - i(AOA(xi)A(x2))e-i(xi,x2) 
-(AOA(xi))e-i(xi,yi)(r(yi)-p(yi)) 

+ i(AOA(xi)A(x2))e-i(xi,yi)e-i(x2,y2)(r(yi)-p(yi))(r(y2)-p(y2)). 
As in the first order case, we now have a look at O = A(z) and O = A(zi)A(z2): 



(66) 



(A(z)PP) 



(PP) 



eq. (66) 



- (r(z)) 



(67) 



which is the same result as in the first-order expansion, since all the additional terms due to the second order vanish; 
next 



(A(zi)A(z2)PP) 



(PP) 



eq. (66) 



:r(zi)-p(zi))(r(z2)-p(z2)). 



(68) 



2nd 



Thus, in the second-order Kamlah expansion, both (A(z)PP)/(PP) and (A(zi)A(z2)PP)/(PP) are reproduced exactly 
(cf. sec. |lVE| ) . This gives rise to the hope that the energy, too, will be projected correctly, since it depends both 
on (r(x))(r( y)) and on (r(x)r(y))c. The detailed derivation of the simplified expression for the projected energy is 



given in app. C 1. The main lessons are that, under the proposed assumptions, the magnetic energy does not obtain a 
correction from the projector (since it depends on transversal degrees of freedom only) whereas the correction of the 
Kamlah expansion subtracts just the longitudinal part of the electric energy; one obtains eventually a concise formula 
for the projected energy: 



(HPP) 



(PP) 



2nd 



= {H) - -(Pj;(y,x)n,(x)n,(y)) + -p(x)GA(x,y)p(y). 



(69) 
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Note how different the outcome is from the first-order calculation: the projected energy is completely independent of 
the longitudinal parts of the electric field (both the expectation value and the correlator), and in turn they cannot be 
determined variationally. In the wave functional, the longitudinal parts are totally undetermined, which is what one 
wants since their value is prescribed by the projector. Thus we can conclude that, at least in electrodynamics, the 
second order Kamlah expansion is mandatory for a consistent formalism. 

One amusing point might be noted at the end of this section; if one restricts the dimensionality of space-time to 
(1-1-1), there are no transversal components of A, 11; thus the projection to second order as carried out by the 
Kamlah expansion leaves behind only the interaction of static charges. Though we have not checked explicitly (by 
replacing Fourier integrals by sums etc), eq. (|6^) seems to suggest strongly that if we consider the theory on a compact 
spatial manifold in (1-1-1) dimensions, we would also keep - in addition to the energy of the charges - the spatially 
constant part of 11 (since this part is annihilated by the derivatives in the longitudinal projector, and is thus not 
subtracted), just as one would expect. 



E. General Observations 



A certain pattern could be observed when going from the first- to the second-order expansion: in the first-order 
expansion the expectation value of a single Gauss law operator between projected states was corrected, but not 
those of more operators (we checked the expectation value of A(zi)A(z2) in eq. (p3|)). In the second-order expansion 
the expectation values of up to two Gauss law operators between projected states were corrected, but not those 

of more operators (it is easy to convince oneself that (A(zi) A(z2)A(z3)P'')/(P'') ^ipj — (A(zi)A(z2)) (r(z3) — 
p(z3)) -|- permutations). Thus, we would expect that in order to have the expectation values of up to three Gauss 
law operators corrected we need a third-order expansion etc, although we have not checked this explicitly. To what 
extent this pattern extends also to other operators is currently not known. The Gauss law operator is in this context 
obviously useful as a first check on the accuracy of the method, since one knows what to expect of the projected 
expectation value beforehand. One can give a further argument in support of the pattern found here: If we want to 
project correctly the product of n Gauss law operators we have to obtain a product of n charge distributions p from 
the Kamlah expression for the projected operator. We have seen that the coefficients Ai are independent of p so the 
only place where p does appear is in the terms ((F) — pY (where for simplicity we haven't written out the product). 
But the highest power i that we obtain is equal to the order of the Kamlah expansion, i.e. a Kamlah expansion to 
second order includes at most terms cx p^. Thus the operator F^ can never be correctly projected to give p^ in the 
second order Kamlah expansion and a third order treatment would be called for. 



V. APPLICATION TO ABELIANIZED YANG-MILLS THEORIES 



A. Difficulties with External Charges 



After the experience with the Kamlah expansion in electrodynamics, we want to try to apply it also to the Yang- 
Mills case. We will see, however, that we cannot incorporate charges as simply as in electrodynamics, and an explicit 
computation does not transcend the leading order in perturbation theory (although a possible background field is 
treated to all orders). 

In order to discuss why we cannot include charges, the first order of the Kamlah expansion is sufficient; we don't 
consider the first order any further than necessary for this purpose, since its shortco ming s were already obvious in 
the discussion of the electrodynamic case. Again we start from the expression (cf. eqs. (ISU |3^)) 



{He-' f '^°^'") = An{e-' f '^°^'") + A?(y)(A''(y)e~* '' '^°^'") 



(70) 



where we have used the abbreviation A°(x) = r"(x) — (F"(x)); later on, we will also use^ 6°'{'x.) = p°(x) — (F°( x)). 
In contrast to the previous section, (5"(x) is now an operator. Using the procedure that is outlined in app. CI, we 



^ Here, a little care has to be taken w.r.t. the placement of gs. If we use 'perturbative scaling', cf. app. IaI the generator of gauge 
transformations in the gluonic sector is not —F", but —ir"^. Thus the projector is in principle J'Dipe~^I^f ^ = jDipe'I^ ~3P ). 
However, this explicit scaling will only be used later on when we have already done away with the charges; we use the 'non-perturbative 
scaling' while we are developing most of the formalism, and re-install explicit factors of g later on. 
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obtain 

(i/POIp^ {H){PP) + (YM|i?r'^(x)|YM)(e-i)'^^(x,y)(<5^(y)P''), (71) 

where we have used that r''(x)P'' = p"-{x)PP, and that our states |) have the form |fermion} |YM) (This is actuaUy 
a restriction - cf. eq. (|lj) - which is of no importance to the argument made here). At this point we see aheady 
that the results of electrodynamics do not generalize straightforwardly: there the fermionic states could be chosen as 
eigenstates of all charge operators p. In Yang-Mills theories this is not possible (apart from states which are either 
non-charged or where the charges form locally colour-singlets), and therefore we cannot write 

(<5^(y)P'') = (<5''(y))(P''). (72) 

Thus, we cannot get rid of all expectation values containing the projector simply by considering the normalized 
projected energy expectation value {HPP) / {PP) . This, however, was the main appeal of the Kamlah expansion, and 
in order to rescue as much as possible of it, we will consider only uncharged states in the following: states that satisfy 
(/|p°(x) = and p''(x)|/) = 0. In this case, the projector is reduced to its purely gluonic part, and therefore called 
P in the remainder of this section. 



B. Kamlah Expansion to Second Order 



We will now turn to the second order Kamlah expansion. While determining the coefficients Aq, Ai, A2, we still 
work with general (possibly charged) states; it will turn out that the coefficients are actually independent of whether 
we have a charged state or not. We have set E = throughout the whole calculation. Then, we start from (cf. 
eqs. (|1|,|2|)) 

{He-^^^^^^} = Ao(e-' f + A?(y)(A'^(y)e-'/ '^°^") -I- (y, z)(A"(y) A''(z)e-'^ / '^°^°). (73) 

It seems to be quite sensible to require that A2^{y,z) should be symmetric under interchange of all indices, i.e. 
A2''{y,z) = A2"'{z,y), since (a) most of A°(y)A''(z) is symmetric apart from r'^(y)r''(z) ^, (b) the antisymmetric 
part of A°(y)A''(z) is — |/''''^r'^(y)(5yz; and should therefore rather be assigned to the term multiplying Ai from a 
systematic point of view. In the course of the calculation, we have to perform two approximations: first, we have 
to restrict ourselves to Gaussian states (as in electrodynamics), cf. eq. (^). However, we have to make a second 
approximation, which is more severe: we will restrict the (gluonic part of the) Gauss law operator to its 0{g^) part. 
In order to facilitate this perturbative treatment, we start with the 'non-perturbative scaling' scheme (cf. table | in 



app. A 5) and then replace A by A -I- ga. Then, a is the new dynamical quantum field with (a) — 0, and jj^ is in 
the following denoted by 11. The pair (a, 11) then follows the 'perturbative scaling' scheme. Note that in this scaling 
scheme we continue to call D,°''(x)n^(x) the (gluonic part of the) Gauss law operator, even though the generator 
of gauge transformations in this scheme is — iDf''(x)nJ(x). The Gauss law operator is now expanded w.r.t. g as 
follows: 

r'»(x) = Df (x)n,^(x) - (Df (x) - g/-\^(x))n^(x) = f '^(x) - .9/-''<(x)n,^(x) 

= f-(x) + 0(g), (74) 

i.e. r°(x) = Df''(x)n^(x), where D°^(x) is the covariant derivative in the background field A, explicitly 'Df''{x) = 
y.gab _ jfacb^.(x) . This then will allow to use factorizability in a fashion identical to the case of electrodynamics. 
However, this approximation has consequences; the first one is that we have reduced the non-Abelian SU(N) symmetry, 
cf. eq. (^, to an Abelian {/(l)*^^ symmetry since: 

[fnx),f''(y)]=0. (75) 

The second one has to do with the admissible background fields A; it is clear that we cannot expect accuracy 
of the energy expectation value to orders higher than 0{g'^)] thus we approximate the Hamiltonian by its 0{g'^) 



6 Note that [r", !•*'] ^ 0, cf. cq. {|). 
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contribution^. The kinetic energy remains unchanged, but the magnetic energy is simphfied quite dramatically if we 
insert the decomposition A°(x) — A,°(x) +gaLf{x) into p-B"(x)B"(x) and neglect all terms of higher order than g'^ : 



-B^(x)B^(x) 

9 

1b^(x)B^(x) + ^B^(x)e,,,fe,(B;^^^(x)5.,Ja^;(zi) 



1b^(x)B^(x) + |j,;nx;zi)a^;(zi) + A4^^^g(x;zi,Z2)a^;(zi)a-(z2). 



(76) 



(77) 



For eqs. ([T^, ^7|) to be valid we do not have to integrate over x. The restriction on A comes about from the 
requirement that the approximated Hamiltonian should commute with the approximated Gauss law operator, otherwise 
the (implicit) assumption of the Kamlah expansion - that in the beginning the Hamiltonian commutes with the 
projector - is incorrect. That one obtains a non-trivial condition at all results from the fact that, to 0{g'^) of the 
commutator of the full quantities, one gets also a cross term from the 0{g^^) part of B^ with the 0{g^) part of F. 
Thus we can express the requirement of gauge invariance of the approximated Hamiltonian in two different ways: 



(a) the term of order 17 ^ in eq. (77) should vanish, since then there can be no cross term with the 0{g^) contribution 
to F. With a partial integration, we see that this condition can be written as 



6,,,,,D;f (x)B?(x) = -D^\^''(x)F;^,Jx) = 0, 



(78) 



which is nothing but the classical equation of motion for the (static) Yang-Mills field, since F is the field strength 
tensor. 

(b) the commutator p-[B^,f] should vanish. We can compute the commutator and obtain 



[B,?(x)B^(x),P(y)] ^ 2Df;'^(y)M,%^;(x;y,Z2)a^^Jz2) ^ 0, 



,fblC2l 



9' 



(79) 



where we have used the fact that M - as defined by eqs. (^, |7^) - is symmetric under interchange of all its 
indices. This is of course equivalent to condition (a), and thus to the requirement of the background field 
fulfilling the classical equations of motion for static fields^ . 



The computation is given in some detail in app. C 2 c. Here we want to mention that in the course of this co mput ation 
a further restriction has to be made of the class of admissible kernels G^^ which is stated explicitly in eq. (|C26| ). We 
then obtain an overall correction to the electric energy 



(80) 



where H^ is the (generalized) longitudinal projector as defined in app. A 6. Thus, the electrical energy neither 

contains D"*e^ nor the (generalized) longitudinal component of G~^. Since by requirement of approximate gauge 
invariance the magnetic energy is transversal as well, this is true for the total energy, and to 0{g'^), we can write 



{HP) 



2nd 



i J d3:,,d^X2(HT)ttHx2,Xi)(n?^nxi)ng(x2)) 



d-^x 



,B«(x)B^(x) + tr(A^G: 



(81) 



This approximation is also sensible from another point of view: the 0(9") contribution of the Hamiltonian is at most quadratic in the 
operators A and 11 and from the case of electrodynamics on e exp ects that the second-order Kamlah expansion treats correctly only 
(projected) expectation values of up to two operators, cf. sec. IV E. 
* In the electrodynamical case there is no 0{g^) contribution to the Gauss law operator and hence the whole discussion is not necessary. 
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with tr(7WG'xy) = ((D(x)D(x) - B(x))°j'G'5^(x, y)), where we have used the abbreviations D^^x) = 

ab 



eifcjD^^(x), (x) = B^(x)/'^'^^eifcj; for IIt cf. eg. ( |A26| ) . Writing the kinetic energy also in terms of G ^ and 
e, the energy is given as 



approx 



{HP) 



(P) 



2nd 



\ I cPx.Sx^ {(nT),tt(^2,xi)<^(xi)e?^^(x2)} + iTr(Gj;^) 
4 / "^'"^ {^B^(x)B^(x) + tr(A^Gxy)|y^x| , (82) 



where is defined by eq. (C26), and the trace Tr furthermore runs over colour, spatial and position indices. We see 
that the projected energy is in fact independent of the (generalized) longitudinal components of e and G^^. This is 
what we would expect from the projector: we treat it to 0{g'^), and we obtain an energy functional that is compatible 
with the energy functional of a state that is annihilated to 0{g'^) by the Gauss law operator, cf. app. Perturbation 
theory thus seems not to invalidate the treatment of the Kamlah expansion. It is interesting to note that eq. ( |8^ ) 
is (at least formally) identical to the improved energy functional obtained in in the case of magnetic background 
fields, to the same order in perturbation theory as considered here. One difference, however, should be noted: in our 
treatment the condition that A has to satisfy the classical equations of motion (and thus has to be transversal 
in the generalized sense) occurred naturally due to the requirement that the approximated Hamiltonian should still 
commute with the approximated Gauss law operator, whereas in |^ it is apparently assumed that is transversal in 
the generalized sense. In the application considered in [^, the Savvidy vacuum, the condition is of course fulfilled. In 
other cases, as e.g. the constant magnetic field stemming from a non-commuting gauge potential p^ , |2^ , this need 
not be the case. 

As a last topic in this section we want to consider the value of the energy if we consider it at the stationary point of 
e and G. We first consider e: 



5 



<5e-(x) 



K7^ro. = J d'x, {(nT)r/(x,Xi)e?^Hxi)} = 0, (83) 



thus the transversal component of e has to be zero at the stationary point. Next we consider Gtt (since Ai is 
transversal, the magnetic energy also contains only Gtt)- For the variation w.r.t. Gtt we need the relation 

Tr(G-i) = -(G4g4)^^7(zi,Z2). (84) 



^(Gtt)^/ (Z1,Z2) 

After a short calculation we end up with 



KllioAM - ^Tr((DD _ B)^) + ^ 1 d3xB«(x)B-(x). (85) 

Thus we see that at the stationary point, the difference to the treatment of [Q, where the projector was ignored, 
lies only in the requirement that here A has to satisfy the classical equations of motion. The important point of the 
Kamlah expansion is, however, that the energy functional depends only on the transversal degrees of freedom even 
away from the stationary point. 



C. Treatment of Charges to Lowest Order in g 



Here we only want to insert an additional point that is inspired by the observation that, in the leading-order 
perturbative treatment, the gauge group SU{N) is reduced to a direct product of U{\) groups. In that case one can 
obviously also require the charged states to be simultaneous eigenstates of all charge operators, and can therefore 
factorize 

{5\y)P'')^5\y){PP), (86) 

where the &^{y) on the RHS is no longer an operator but simply a c-number function. At this point, we should make 
clear that in the 'perturbative scaling', cf. app. |^, that has been used here, 5 reads 



5''(x)=gp"(x)-(f"(x)), 



(87) 
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since Gauss' law implies Q""]) = (— + = — > r°|) = gp°'\)- For this kind of quasi-Abelian charges, we can 
immed iatel y obtain the projected energy to second order in the Kamlah expansion by replacing (F") by (F°) — gp" 



in eq. ( |C2l|) : 




= (H) - -(Ai7A=Kxi)A=^(x2))e^^^Hx2,xi) 



-(Ai?A^nxi))e^^^nx2,xi)((r^(x2)) - .gp=^(x2)) (88) 
+ i(AiJA^i(xi)A=^(x2))e^i'^nxi,yi)e^^''Hx2,y2) X 



We see that it differs from the energy projected onto the chargeless sector by 



-g{ {{AHA^- (xi) A=^ (X2))e=^'^^ (xi , yOG^^"^^ (x2, y2)(f (yi))) 
-(Ai^A^i(xi)}e^^'^Hxi,y2)}p'^Hy2) 
+ ^ (AiJ A^^ (xi) A^^ (X2))e^^'^^ (xi , yi)e^^'- (X2, y2)/^ (yi)p''^ (ya). 



(89) 



The first term can be easily calculated to give zero ^, thus there is no term linear in p that contributes to the energy, 
whereas the second term simplifies to 



We see that in leading order perturbation theory, the charges interact via a potential given by the Green's function of 
the background field-covariant Laplacian. Three comments are in order: first, since the Green's function depends on 
the background field, the variational equations for A are not independent of the charge distribution. This is different 
from the electrodynamical case where radiation (transversal terms) and static charges were completely decoupled. 
Second, we see that obtaining the (3 function from corrections to the Coulomb potential (which would result for 
A = 0) is on this level of approximation not possible; it seems almost surprising that we obtained the correct leading 
order potential, since our projection scheme only attempts to give the correct energy to 0{g'^). For the first quantum 
corrections to the interquark potential, we would need the energy up to 0{g'^). Third, the fact that the inter-charge 
potential is given by the Green's function of the covariant Laplacian strongly reminds one of the result of [ p^ where 
a rigid rotator model was considered which was also inspired by ideas from nuclear physics. 



In this paper we considered the Kamlah expansion, which is an expansion of the expectation value of the Hamil- 
tonian between projected states in powers of the symmetry generator. This is a technique successfully employed in 
nuclear physics, but it is not clear a priori that the same can be done in a field theory with a local symmetry since the 
applicability in nuclear physics is based on the fact that the states are strongly deformed, which comes about usually 
from many particles participating in that deformed state. 

We tried the Kamlah expansion first for electrodynamics with external (static) charges, where we found that the first 
order of the Kamlah expansion is inconsistent with its assumptions, since at the stationary point we find a state 
that is not deformed. Going to second order, however, we found exactly the properties that we expect of a projected 
energy functional, namely that it is independent of the longitudinal parts of the (11 11) correlation function and the 
longitudinal part of the expectation value of 11, and the correct Coulomb energy appears in the projected energy. 



We have already indicated above that the two apparently very different moments of inertia Xsc,Iy are identical in a perturbative 
treatment to leading order in g; this is the point where one observes this most clearly. For the evaluation of the expressions used here, 

we have used a different form of the longitudinal projector, (11^;,)^^^'^^-'^ (x2, xi) = ^D^^*^ (x2)D^^''^ (xi)G^''^ (xi , X2)^ , where again we 

don't integrate over xi,X2. 





(90) 



VI. SUMMARY AND CONCLUSIONS 
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This result motivated us to consider Yang-Mills theories, too, in this framework but we restricted ourselves to the 
leading order in a perturbative expansion in powers of the coupling constant g. Since we have approximated both 
the Gauss law operator and the Hamiltonian, we obtain consistency conditions from the requirement that the ap- 
proximated Hamiltonian still commutes with the approximated Gauss law operator. Effectively, this restricts the 
choice of possible background fields to those which satisfy the (static) classical equations of motion. This condition 
makes the operator that multiplies the part of the potential energy quadratic in A transversal w.r.t. the covariant 
derivative in the background field. Since the projection subtracts off the longitudinal part (in the generalized sense) 
of the electrical energy, we end up with an energy functional that is independent of the longitudinal parts (in the 
generalized sense) of the parameters of the Gaussian wave functional. At the stationary point of all parameters save 
the background vector potential, the energy functional looks precisely like the one-loop functional that one obtains 
in the mean field considerations without taking into account the projector. Of course, this is so only formally, since 
on the one hand the consistency condition alluded to above ensures that the energy depends only on the transversal 
degrees of freedom (thereby reducing three polarization states to two), and, on the other hand, the longitudinal parts 
of the parameters in the Gaussian wave functional are undetermined - instead of being set to zero - just as one would 
expect from a projected energy functional. We then also considered how external charges can be included into the 
Kamlah expansion. Up to now, we have been able to allow for them only to leading order in perturbation theory, 
since only then the charges can be treated in a quasi- Abelian manner. 

Three investigations seem worthwhile: first, one has to investigate whether the Kamlah expansion is a truly non- 
perturbative expansion or whether at some stage an implicit expansion in powers of g takes place. In this context, one 
should also investigate further the 'pattern' found in our discussion on electrodynamics, namely that one seems to need 
the Kamlah expansion to n*'' order to project correctly terms containing n operators. If this would be true in general, 
one would need the fourth-order Kamlah expansion for a correct projection of the full Yang-Mills Hamiltonian. 
Second, one should try to find out whether the Kamlah expansion can be carried out also to higher orders in pertur- 
bation theory since this would lift the restriction to background fields that satisfy the classical equations of motion. 
This should be considered also in the presence of external charges, which up to now can only be dealt with classically. 
If that can be done successfully, one should face the problem of computing the one-loop contribution to the inter- 
quark potential. If this reproduces the standard results, one could also try in a third investigation a non-perturbative 
evaluation of the second-order Kamlah expansion (provided the Kamlah expansion is truly non-perturbative and a 
second-order treatment is sufficient), although this will be a quite difficult task, since in our evaluation we relied 
strongly on the factorization properties; these, however, cannot be applied so easily in a non-perturbative framework, 
since the full Gauss law operator contains products of two operators at the same point in space. 
Thus, renormalization problems, which we have ignored completely in this paper, will have to be taken into account 
in a non-perturbative manner before one can perform such a non-perturbative evaluation of the projected energy 
functional. 

Of course, when considering such non-perturbative extensions one should keep in mind that the restriction to Gaussian 
states has tied us already to the Hartree-Fock approximation, which is basically the two-loop approximation of the 
total quantum energy. 
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APPENDIX A: CONVENTIONS 
1. Indices and Summation Conventions 

In this paper, four different kinds of indices appear: colour indices, spatial indices, position indices and super- 
indices, combining several of the former three kinds as will be explicitly stated wherever they are introduced. To 
each kind of index a certain range of letters is attached: colour indices: a, b,c,.. ., spatial indices: k,. . ., position 
indices: x, y . . ., super-indices: i,j,k,... . 

Unless explicitly stated otherwise, Einstein's summation convention will be used throughout. Discrete indices will 
be summed over, continuous indices will be integrated over. We will also use a uniform notation for S functions, 
regardless of the nature of indices (discrete or continuous); thus, (5xy in this paper corresponds to (5(^)(x — y) in the 
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usual notation for Dirac's S, whereas e.g. S is simply Kronecker's S. 

2. Minkowski Space 

We use as metric for the Minkowski space 

= diag(l, -1,-1,-1), (Al) 

and thus have as four- vectors a'' and — g^.^a'^ 

a^ = (ao,a) , a^ = (ao,-a) (A2) 

with the scalar product 

a ■ b = a^h^ = aobo — a.b. (A3) 

Only for the four-derivative, signs are distributed differently: 

9'' = (9o,-V) , d,^{do,V) (A4) 

with Vi — Sometimes we need a gradient w.r.t. another variable. This we denote by attaching the variable name 
toV, e.g. Vf = 

3. Group Theory Conventions 

The (hermitian) generators of the SU(N) Lie algebra satisfy the commutation relations 

[A'^,A''] =ir''"A". (A5) 
The generators in the fundamental representation(s) are normalized to (1/2): 

tr(AfVdALd) = l^^'- (A6) 
The structure constants satisfy the Jacobi identity 

jabr jcdr _|_ jacr jdbr _|_ jadr jbcr g (A7) 

and are normalized as 

j^abcjdbc ^ j^gad_ (-^g-) 

4. Canonical Quantization 

Throughout this paper we work in the Hamiltonian framework. We use the Weyl gauge = 0; but do not fix the 
gauge any further. The usage of the Weyl gauge implies that the physical states of the system have to satisfy the 
Gauss law constraint. In the absence of external charges (which we assume throughout this appendix), the condition 
reads 

r"(x)|V')=0 (A9) 
where r°(x) = Df (x)n^(x) with Df (x) = Vi^"'' - r"''A^(x). (AlO) 

D°''(x) is the covariant derivative in the adjoint representation. If we consider a system with external charges, the 
Gauss law operator obtains additional contributions that stem from the external charges, as is outlined in sec. ||. 
In order to quantize the system, we do not have to resolve Gauss' law, since even a system with redundant coordinates 
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may be quantized We work in the coordinate representation, i.e. A is realized multiplicatively, whereas the 

canonical momenta are realized as derivative operators: 

(A|n- (x)l^) = i^^(AlV^). (All) 

This choice ensures the satisfaction of the equal-time commutation relations automatically 

[A,«(x),n5(y)]=*rt,<5xy. (A12) 

Factors of g appear in seve ral q uantities at different places according to the different conventions possible as will be 
discussed at length in app. A 5. However, contrary to some conventions used in the literature, the meaning of A, 11 
is always the same irrespective of wh ether we choose 'perturbative' or 'non-perturbative scaling', especially no factors 
of g enter the commutator eq. (A12). 

The wave functionals we deal with in this paper are given in the A-representation 

il)[A] = (A|V'). (A13) 

In the language of wave functionals the fact that the choice of Weyl gauge implies that has to be annihilated by 
the Gauss law operator is translated as 

%Ij[A^] = ijj[A] (A14) 

for topologically trivial but otherwise arbitrary gauge transformations U . In this context, A*^ denotes the gauge 
transformed field 

(A,"A")^(x) = [/(x)A^(x)A"[/^(x) - ^^7(x)V,C/t(x) (A15) 
with J7(x) = e*"^°^^^^" and A" are the generators of SU(N) in some representation. 

For further details on the canonical treatment of Yang-Mills theories, cf. [Q; more details on the Schrodinger picture 
can be found e.g. in (2^ G8|. 



5. Factors of g 

In Yang-Mills theories one has basically two options concerning where one wants to put the coupling constant, 
either in front of the action (J non-perturbative scaling'), or in front of the commutator term in the field strength 
{'perturbative scaling'). In table | we give a short list concerning which convention leads to which placing of factors 
of g in other quantities of interest. 





non-perturbative scaling 


perturbative scaling 


covariant derivative 


-Dp = - iA^ 


= df, - igA^ 


field strength 


Fi_iiy = Oi_lAi^ Oi/A/^i Z[ylp, Ajy] 

= i[D^,D^] 


F^„ = d^A^ - d„A^ - ig[A^,, A^] 


action 


S^~^F^^,{x)F^^^{x) 


S = -\F^^,{x)F^^-'{x) 


electrical field 


T^a 77a 


E? = FS, 


magnetic field 


= (V X A)? ir''"(A' X A")i 


B? = -\e,jkF''^'' 
= (V X A)t - f /"''"(A' X A'^); 


momenta 11'' = Sttt 


-pra 1 T-^a 1 -rpa 


n? = i^, = -E? 


Hamiltonian 


H = ^nnx)n?(x) + ^B?(x)B?(x) 
= 2|7(E?(x)EKx) + B?(x)BKx)) 


H = in?(x)n?(x) + iBnx)B?(x) 
= i(E?(x)E?(x)+B?(x)Bnx)) 


wave functional of 
'free' theory 
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non-perturbative scaUng 


perturbative scaling 


gauge transformations 






generators 




[i(-r«),i(-r^)] = 


finite gauge trafos 
(gluonic part) 







TABLE I: Placement of the coupling constant g in the 'non-perturbative' 
and the 'perturbative' scaling scheme. We have used the shorter form 
n for r"(x). 



6. Generalized Projectors 



In this appendix we want to give some notions on the generalized projectors as they are used in this paper. First let 
us note that the projection operators are thought of as bilocal objects, i.e. they depend on two spatial coordinates. 
Consider e.g. the (generalized) longitudinal projector (11^): 

(ni)«''(x,y). (A16) 

As is to be expected for a projector, (11^,) is idempotent (if we include an integration over the double continuous index 
as well): 

(ni)^''(x,y) = J d^z(ni)^,?(x,z)(ni)^.5(z,y). (A17) 
The (generalized) longitudinal projector (11^) is a symmetric operator in the sense that 

J d^xd^yfn^)iIiL)t^{^,y)9',{y) = J d^xd^ygti^){IlL)t^i^,y)f^iy), (A18) 

where we have dropped boundary terms from two partial integrations. The latter property implies that it can be 
diagonalized, and from eq. (A17) one can conclude that its eigenvalues are 0, 1. The generalized longitudinal projector 
can be given in different forms; for this it might be useful to recall that the ordinary longitudinal projector also 
can be given in different forms that are equivalent under the usual circumstances (i.e. all functions of y are found to 
the right 10 of (P^)„ (x,y)) 

(P^).,(x,y) = -fvfGA(x,y))vy (A19) 



VjvrGA(x,y)j (A20) 
rf'fc^^.k.(x-y)ll^^ (A21) 



(27r)3 k2 

where G'a(x, y) denotes the Green's function of the Laplace operator, — ViViGA(x,y) = ^xy. Under the given 
circumstances all three definitions are equivalent; however, only the first one really has to 'act to the right' in order 
to make sense. The (ordinary) transversal projector can also easily be defined as 

(P^)„(x,y) = 5,,6^y - (P^)y (x,y). (A22) 



In other words, for expressions of the type /i(x)(P^ )ij (x, y)gj (x, y) (with a sum or integral over all double indices) all these definitions 
are equivalent, otherwise they will usually be not. 
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For the generalized longitudinal projectors we can now give corresponding definitions: 

(ni)^^(x,y) = -(Dr^(x)Gl^''nx,y))D^'^'(y) (A23) 
= (Df^(y)Dr(x)G2,^''^(x,y)) (A24) 

n 

where G''^'"{x,y) is the Green's function of the covariant Laplacian (DD)=i'^i (x), i.e. -(DD)'=i'^i (x)G^i''^ (x, y) = 
(5^i''^(5xy, (^„)f(x) are the normalized eigenf unctions of (ni,)^j'(x, y) with eigenvalue 1, and we do not integrate over 
X, y. We assume that the covariant Laplacian has no eigenvalues zero. Here the corresponding generalized transversal 
projectors are given by 

(nT)^^(x, y) = S^.S^'S^y (nL)^^(x, y). (A26) 

APPENDIX B: RESULTS FROM PERTURBATION THEORY 

In this appendix we want to address the following problem: assume that we start with a wave functional that has 
Gaussian form. Since this wave functional is not annihilated by the Gauss law operator r°(x), we multiply it by a 
polynomial. How do we have to fix the coefficients of the polynomial s.t. Gauss' law is satisfied to a given order in 
the coupling constant g ||29|] ? 

In order to proceed, we use the same procedure as in the main text: start out from 'non-perturbative scaling', and 
then replace A by A + f;a so that the background field does not decouple in the perturbative limit. The covariant 
derivative can then be split into a part containing the background field (which is completely of 0{g'^)) and a part 
that contains the fluctuation part and an explicit factor of g: 

Df (x) = Df (x) - g/'^"\"(x) with Df (x) = \/,S''^ - /'^"^A^(x). (Bl) 
Thus the Gauss law operator reads (cf. the discussion in sec. |VB| ) 

r-(x) = Df (x)i-i- = iDf (x)-i- -5r=^<(x)i-i-. (B2) 
I oa°(x) I oai(x) » <ja°(x) 

The wave functional is of the form 

where we have used a (hopefully) obvious shorthand notation together with the abbreviation x — ^G^^ — 2iS, and 
P(a) indicates a power series in g: 

oc 

P(a) = E5"p(")(a). (B4) 

If we now require 

r'^(x)* = + 0(5^+1) (B5) 
we obtain a recursion relation between the different P'^"^s: 

N N 

^ ^ AP(") ^ /iP("-i) 

+ J2 — ^ - J2 f'^r'^i^ = 0, (B6) 
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where it is sufBcient to terminate the sums at n = A'" since we require Gauss' law only up to 0{g^). Note that the 
summation indices start at different lower values. Thus we can study the 0(5°) and 0{g^) cases separately^^: 

0(.g°): Df((-xaX'+ze^)-0. (B7) 

Since this equation has to be valid for all a, one actually obtains three conditions since D does not mix real and 
imaginary parts: 

Dx = ^ Df (x)(G-i)^/(x,y) - and Dr(x)S^/(x, y) = (B8) 
De = 0, i.e.Df (x)e^(x) = 0. (B9) 
Note the similarity with the QED condition of having a purely transversal kernel. 

Carrying this procedure further to higher orders in g will lead to certain difficulties, e.g. that ip will be annihilated 
only if one takes a to stem from a subspace satisfying the background field condition Da — 0, thereby invalidating 
the whole concept, since our states should be such that Gauss' law is satisfied for all a without further gauge fixing. 
However, to 0(5"), the concept works, and is compatible with the results from the Kamlah expansion as seen in sec. 

The main lesson from this section is that a pure Gaussian state can satisfy Gauss' law only to order g". This is the 
reason why we can formulate an exactly gauge invariant Gaussian ground state for electrodynamics. 



APPENDIX C: DERIVATIONS 

This appendix contains some of the technical details of the different derivations. 



(CI) 



1. Electrodynamics, Expansion to Second Order 

a. Generalities 
Here one starts from the expression (cf. eqs. (^ ^2|)) 

+A2(y,z)({r(y) - (r(y))}{r(z) - (r(z))} e'/^(r-p)). 
The coefficients will be determined as follows: one needs three equations, which are obtained by 

(a) setting (f> = 0, 

(b) performing a functional derivative w.r.t. 4'i'^i) a-nd setting = afterwards, 

(c) performing two functional derivatives, w.r.t. (/)(xi), 0(x2) and setting </> = afterwards. 
After a bit of manipulation, one obtains 

Ao = (ff)-A2(y,z)(A(y)A(z)), (C2) 
Ai(y) = ((FA(xi))-A2(yi,zi)(A(yi)A(zi)A(xi)))e-i(xi,y) (C3) 

and for A2 one obtains an even less pleasant equation 

(ifA(xi)A(x2)) - (ff)(A(xi)A(x2)) - (ifA(z))e-i(z,y)(A(y)A(xi)A(x2)) 
= A2(y,z)[(A(y)A(z)A(xi)A(x2)) - (A(y)A(z))(A(xi)A(x2)) 

-(A(y)A(z)A(yi))e-i(yi,y2)(A(y2)A(xi)A(x2))l, (C4) 



11 We set P(0) 

= 1, since in the end it will be fixed by the overall normalization anyway; it has to be nonzero, s.t. the state does not 
vanish as <; — ► 0. 
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with A(x) = r(x) — (r(x)) and 8 -'^ is defined via Q "'^(x, y)(A(y)A(z)) = (5xz- Eq. (C4) simplifies significantly if one 
exploits sonic of the factorization features of Gaussian states, 

(A(x)A(y)A(z)) = (C5) 
(A(y)A(z)A(xi)A(x2)) - (A(y)A(z))(A(xi)A(x2)) + (A(xi)A(z))(A(y)A(x2)) 

+ (A(x2)A(z)}(A(xi)A(y)). (C6) 

This allows to rewrite eq. ( |C4| ) as 

(i/A(xi)A(x2)) - (i/)(A(xi)A(x2)) 
= A2(y,z)((A(y)A(xi))(A(z)A(x2)) + (A(y)A(x2))(A(z)A(xi))). (C7) 

In order to obtain A2 explicitly, one can again use 0~^: 

((A(y)A(xi))(A(z)A(x2)) + (A(y)A(x2))(A(z)A(xi)))ie-i(xi,yi)e-i(x2,y2) 



), (C8) 



_ 1 
2 

giving (since A2(yi,y2) = A2(y2,yi)) 

A2(yi,y2) = ((i?A(xi)A(x2)) - (i/)(A(xi)A(x2))) ie-i(xi,yi)e-i(x2,y2). (C9) 
This result leads directly to eq. ( |65|) in the main text. 

b. Explicit Form of Energy Functional 

We will now turn to the evaluation of the energy. In addition to the terms we have computed in the context of 
the first order Kamlah expansion, we need to compute (Bi(x)Bj(y)r(zi)r(z2)). If we again - as in the first order 
calculation - take either E = or G^^ and E translationally invariant, the matrix element factorizes, 

(B,(x)B,(y)r(zi)r(z2)) = (B,(x)B,(y))(r(zi)r(z2)) (CIO) 

which is again plausible, since there should be no corrections to the magnetic part of the energy as it depends on 
the transversal degrees of freedom only. The next ingredient we have to compute for (Ai? A(xi)A(x2)) is the electric 
energy, which gives a rather simple contribution, too. We can evaluate it using the same factorization properties that 



underly eq. (C6), and obtain altogether: 

(AFA(xi)A(x2)) = (r(x2)n,(x)},(r(xi)n,(x)},. (cii) 



-(Ai/A(xi)A(x2))e-i(xi,x2) - -(Pf,.(xi,x2)n,(x2)n,(xi))„ (C12) 



This expression enters the correction of the mean-field energy 

i(Ai/A(xi)A(x2))e-i(xi,x2) = i 

where P^(x2,Xi) denotes the longitudinal projector as defined e.g. in eq. ( [A19| ); translational invariance, that has 
been assumed before in the first order calculation, implies that the electrical energy can be written as 

i(n,(x)n,(x)) = i(P^,(y,x)n,(x)n,(y)) + i(Pj:(y,x)n,(x)n,(y)) (C13) 

and thus the second order Kamlah expansion precisely subtracts off the longitudinal (gauge variant) part of the kernel. 
The other term needed to compute (iJP'')/(P'')|2„is even simpler, 

i(Ai7A(xi)A(x2))e-i(xi,yi)e-i(x2,y2) = GA(yi,y2). (C14) 

Similar to the first-order discussion, we subtract off the part (nf )(nf ) (thus completing {P^.^U,Uj)c to {Pf-U^Uj)) 
and add it back in afterwards. Thereby, we obtain eq. (|69|). 



In this context there is no problem with the other definitions eqs. (A20, A21) 
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2. Yang-Mills Theories, Expansion to Second Order 

a. Generalities 



We follow the same procedure as it was outlined in the case of electrodynamics in app.|C_l| From eq. (|7^) we obtain 
by this procedure again three equations: 



{H) 



(i/-{a^^(xi),a^nX2)}) 



Ao+Af(y,z)(A'^(y)A^(z)), 
Ao(t/^(x))+A?(y)(A'^(y)e^(x)) 

+Af(y,z)(A''(y)A''(z)g^(x)) 

Ao(i{a^^(xi),g^^(x2)})+A?(y)(A'^(y)i{g^nxi),^^nx2)}) 
+Af(y,z)(A'^(y)A^(z)i{e^Hxi),a^Hx2)}) 



(C15) 
(C16) 

(C17) 



where in the last equation, eq. (C17), we have used the abbreviation {C/'^^ (xi), (y^^ (X2)} = (y^^ (xi)5°^ (X2) + 
G'^^ {^2)G'^^ (xi)- We then can - again - express Aq, Ax in terms of A2 and unprojected expectation values; the equation 
for Ai is very similar to eq. (C4) in the electrodynamic calculation, but here it is even worse, since [A°(x), A^(y]l ^ 0. 
In the case of electrodynamics, significant progress could be made using the factorization features eqs. (C5, C6). In 
order to mimic the treatment presented there , and thus to make further progress, we have to resort to the perturbative 
expansion of eq. (fsh as outlined in sec. VB, 



h. Factorizations 



For the quasi- Abelian form F" we again have the factorization properties of Gaussian states similar to eqs. (C5, 
1): 



(A''(x)A''(y)A^(z)) = 0, 

(A''(y)A''(z)A^nxi)A^^(x2)) = (A'^(y)A''(z))(A'=Hxi)A^=(x2)) 

+ (A^Hxi)A''(z))(A'^(y)A^=(x2)) 

+ (A^=^(x2)A''(z))(A=^(xi)A''(y)), 



(C18) 



(C19) 



where we have introduced the notation A"(x) = F'^(x) — (F''(x)}. 



Derivation 



We now use the factorization properties to compute A2 (or, if we hadn't required A2 to be symmetric from the 
outset, the symmetric part of A2; but since the Fs commute anyway, the distinction does not play a role): 



4^'nyi,y2) = i[(i?A^Kxi)A=Hx2))-(if)(A=Uxi)A=Hx2))] X 

(e-i)^i'^nxi,yi)(e-i)^^'^^(x2,y2). 



(C20) 



With this, we can give now explicitly the projected, normalized energy functional to second order in the Kamlah 
expansion (in absence of external charges): 



(P) 



(H) - i(Ai/A=^(xi)A^^(x2))(e-i)^^^nx2,Xi) 
-(Ai/A^nxi))(e-i)=i=Hxi,x2)(rMx2)) 

+ i(AifA=nxi)A=^(x2))(e-i)^^'^^(xi,yi)(e-i)^^'^Hx2,y2)(f'^Uyi))(f''Hy2))- 



(C21) 



As in the case of electrodynamics, we may here also insert F into the Kamlah expansion expression instead of H . 
This is due to the fact that f° is gauge-invariant w.r.t. gauge transformations carried out by e*-'^ ^ . If we do this, 
the terms analogous to the second and the fourth term in eq. (C21) (with H replaced by F") are zero, and the terms 
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analogous to the first and third term cancel, so that we obtain the correct result (r°P)|2„d = 0- The computation 
can also be carried out for (f "r''P)/(P) j^nd, and one obtains also the correct result, namely 0; the term analogous to 
the second term gives — (A°A''), whereas the terms analogous to the third and fourth term add up to give — (f '^)(f 
thus cancelling (together with the term analogous to the second) the term analogous to the first term. 
In order to be able to compute the projected energy functional to 0{g^), we now need the following ingredients: 

(A(n,"(x)n?(x))A=Hxi)) = 2(n^(x)A'=nxi))(n,^(x)), (C22) 

(A(n^(x)n,nx))A^Hxi)A^=(x2)) = 2(n^(x)A^nxi))(n,''(x)A=Hx2)), (C23) 
(A(B^(x)B,"(x))A^nxi)) = 0, (C24) 

(B^(x)B«(x)f-(xi)f-(x2)) = ld'xI)'}f^{^,)t)r;'H^2)A'4t{^;K^,^,) 

= 0, (C25) 

where we have set E to zero here and in the remaining computations of this appendix. Also the expression 
/ (i'^a;D^^^''^(xi)D^^^''^(x2)Af;'^^;^^(x;xi,X2) requires a bit of special attention: there is no integration over Xi,X2 al- 
though they appear more than once, but we have to integrate over x, althou gh it seems to appear only once (but this 
is an artifact of the notation introduced in eq. ([ttI)). The last identity eq. ( |C25| ) is true since we had to restrict the 
allowed background fields to fulfil the classical equations as discussed above. In addition, we need a simple expression 
for Q^^; in the electrodynamics calculation we saw that we had to assume translational invariance of G, S. Here we'll 
need something similar, namely we require 

(G-i)^''(x,y) = (n,.)-Hx,xi)(G-i),^^^r.^(xi,X2)(n,.)^^^j(x2,y) 

+(nT)-nx,xi)(G-i)^^^l=^Hxi,x2)(nT)g)(x2,y) 

= (Gzi)^''(x,y) + (G^^)°,''(x,y) (C26) 

with the (generalized) longitudinal projector^'^ 

(ni)^''(x, y) = -Dr(x)Gl^(x, y)Df (y) , where - (DD)-(x)Gl^(x, y) = S'^'S^y, (C27) 

(no integration over x, y) and the (generalized) transversal projector (nT)"j'(x, y) = S'^^SijSxy ~ (nL)fj' (x, y). Thus, 
we don't want to allow for non-zero cross terms IliG^^IlT, H^G"^!!/,. For the explicit construction of 8^^ this has 
the advantage that one can invert the longitudinal and transversal parts individually: 

{Gll)'^^{^,y){GLL)%{y,z) = {Ui^mK,z) (C28) 
and (G4)«''(x,y)(GTT)5^fe(y,z) = (nT)rfe(x,z). (C29) 

We now want to relate the inverse of DiDj(G~^) to DiDjG: we start with the observation that, since G^]^ = 
IIlG~^IIl, we can write (no integration over x, z) 

(Gzi)^''(x,z) ^ Br (x)Df=(z)(iJ,.)=^==(x,z) (C30) 

and correspondingly (again no integration over x, z) 

(Gii)^''(x,z) = Dr (x)Df^(z)(ifi)^^^^(x,z), (C31) 



where Hl,Kl are auxiliary functions that have the important characteristic that they don't carry spatial indices. 



From the requirement / d^y (G^^ )fj'(x, y)(GLL)5^(y, z) = (nL)^^?(x,z) we obtain 



G^^^nx,y) - J d^z ({fit)f^^^{z){HLr'^{^,z)) {KL)''^^Hz,y). (C32) 



Some properties of generalized projectors are given in appendix 
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If on the other hand we start from the defining equation of ^ 



d'y ( ^DfHx)D-^(y)(G-^):f^(x,y) ) (e-^)''=(y, z) = ^-<5: 



-l\bib2 I 
)ij 



we end up with 



4G^^^^(x,y)- J d'z (^it)f)f^'^iz){HLr'^i^,z)) {e-' f^'-iz,y). 



From the comparison of eq. ( C34 ) with eq. (C32), we obtain 

{e-y^^^{z,y)=4{KL)''^^^{z,y) 

and thus conclude 



Dfnx)D™^(y)(e-^)''^''nx,y) = 4(Gii)-(x,y 



be/ 



(C33) 



(C34) 



(C35) 



(C36) 



where again we don't integrate over x, y. We can now compute the correction terms to the Kamlah expansion very 
easily, using the fact that e.g. (Hf (x)A''(y)} = (y)(G£i)f/(x, y): 

1. (A(nf(x)n,"(x))A^i(xi)A=2(x2))(e-i)^i=^(xi,x2) 
= 2(n,.)gr;(x2,xi)(n,^,nxi)ng(x2)). 

2. (nf (x)nf (x)A^i (xi))(e-i)=i^^ (xi, X2)(f (x2)) 
= 2(ni)gr; (X2, xi)(n,^^^ (xi))(n^: (x2)) 

3. (A(nf(x)n^(x))A^i(xi)A'=^(x2))(e-i)^i'^i(xi,yi)(e-i)^^'^^(x2,y2)(f''i(yi))(f'^^(y2)) 

= 2(ni)^^^r;(x2,xi)(n,^^^(xi))(n^:(x2)). 



Inserting these results into eq. (C21), we obtain the results given by eqs. (pO[-|82|) in the main text. 
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